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Abstract 

We discuss a model which gives rise to cosmic self-acceleration due to 
modified gravity. Improvements introduced by this approach are the following: 
In the coordinate system commonly used, the metric does not grow in the 
bulk, and no negative mass states are expected to appear. The spectrum of 
small perturbations contains a localized massless tensor mode, but does not 
admit dangerous localized massive gravitons. All the massive spin-2 modes 
are continuum states. The action of the model, which is an extension of DGP, 
allows to relax the previously known constraint on the bulk fundamental scale 
of gravity. The latter can take any value below the 4D Planck mass. 



1 Introduction 



The DGP model pQ admits the self-accelerated solution [2] which could be used to 
describe [3] the de Sitter (dS) like expansion of the present-day Universe. 

The model exhibits strongly-coupled behavior already at the classical level [I]. 
As a result, the question of stability of the self-accelerated background w.r.t. small 
perturbations, [5] [8], becomes difficult [91 [TU]. 

Luckily, certain non-linear solutions have been found. This is a case for a 
Schwarzschild-like solution, for which the 4D metric was exactly obtained [11], and 
for the Domain Wall solution for which the full 5D metric was found in Ref. [12J. In 
both cases, the mass (tension) of the solution gets screened by gravitational effects, 
and these sources on the self-accelerated background look as if they had a negative 
net 5D mass (tension)0. This suggests that the self-accelerated background may not 
be problem-free in the full non-linear theory. 

What is a root-cause of this behavior? The self-accelerated solution exists only 
for a certain choice of the sign of the extrinsic curvature, and this choice is such that 
it requires a growing metric in the bulk. For instance, in a simplest spatially-flat 
case, the full 5D metric of the self-accelerated solution takes the form [2]: 



where y is the 5th coordinate and H denotes the dS expansion rate of the 4D 
worldvolume (the latter is labeled here by Cartesian coordinates (t,x), and we use 
the 5D coordinate system in which the brane is located at y — 0.). The unusual 
feature of the above metric is that it grows in the bulk, even though the worldvolume 
metric is that of dS space. A linearly growing metric, similar to (CQ), would have 
been produced by a negative tension 3-bran^E The growing metric (0Q) imprints its 
"negative" effects on the brane worldvolume through the extrinsic curvature, giving 
rise to the solutions mentioned in the previous paragraph. 

The aim of the present work is to modify the DGP equations in such a way 
that the new system admits a background that is equivalent to the self-accelerated 
solution on the 4D brane, however, differs from it in the bulk. 

The new solution that we will discuss takes the form: 



In order to obtain this solution one needs to flip the sign in front of the extrinsic 
curvature term in one of the DGP equations, keeping the rest of the equations intact. 
In section 4 we will discuss how such equations can be obtained by modifying the 
DGP action. The latter will have an additional benefit: the constraint [1] on the 
bulk fundamental scale of gravity will be relaxed. This scale, in the present context, 
can take any value below the 4D Planck mass. 

1 In contrast with this, screening of the similar sources on the conventional branch of DGP, 
leaves them with positive 5D mass (tension) [TTJ [TJ] . 

2 This could be in, e.g., 5D Minkowski or Anti de Sitter (AdS) space-time. 



ds 2 = (1 + H\y\f {-dt 2 + e 2Ht dx 2 } + dy 2 , 



(1) 



ds 2 = (1 - H\y\f {-dt 2 + e 2Ht dx 2 } + dy 2 . 



(2) 
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The solution (JSJ) is formally identical to that for a 3-brane endowed with a positive 
4D cosmological constant (brane tension) which is embedded in 5D empty space in 
5D GR [13J. However, unlike the latter, the worldvolume expansion in the present 
case ([2]) is due to modified gravity, while the 4D cosmological constant is set to 
zero. This difference is what is responsible for the modified Friedmann equation, 
and distinct cosmological evolution on the self-accelerated background. 

The spectrum of small perturbations on (T5]) will differ from that of the 3-brane 
with tension, as well as from the spectrum on ([T]). As we will show, the spectrum 
contains one massless spin-2 localized mode, the coupling of which on the brane 
is ~ 1/Mp b and a continuum of massive spin-2 modes with m 2 > 9H 2 /4, with 
suppressed couplings on the brane, due to the 4D Einstein-Hilbert (EH) term. 

The bulk space in both ([T]) and §2§ is locally equivalent to 5D Minkowski space. 
In the chosen coordinate system the solution (TSJ) encounters the Rindler horizon 
at \y\ = H" 1 . However, an analytic continuation beyond this point can be per- 
formed by employing new coordinates (see, e.g., [E]). In that coordinate system 
the brane (with closed spatial sections) can be regarded as a 4D dS bubble that is 
first contracting and then re-expanding in 5D Minkowski space. 



2 Flipped equation 

We briefly review below the DGP equations and self-accelerated solution (QQ), and 
show how modifying one of those equations can give rise to the new solution (j2J). 
Consider the 4D DGP equation written at y = + (we use conventions of [14]): 

G^v - rnJyK^ - g^K) = S-kGnT^x) . (3) 

Here, G^ v = R^ v — g^ u R/2, is the 4D Einstein tensor for the metric ^(x,?/); 
K = g^K^, is the trace of the extrinsic curvature tensor 

K, v = ^ (d y g^ - - V„Ay , (4) 

and V M is a 4D covariant derivative w.r.t. the metric g^ v (x,y). We introduced the 
lapse scalar field N, and the shift vector field [15J: 

9fl5 = iV M = g tlv N v , £55 = N 2 + g, v N» N v . (5) 

Furthermore, the {^v} equation in the bulk, and the {/i5} and {55} equations read 
respectively 

G^=0 for y^O, (6) 
V^ = V,K, (7) 
R = K 2 - K, U K^ . (8) 
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Here G\J denotes the 5D Einstein tensor for the 5D metric qab^x, y) (A,B=0,1,2,3,5), 
and g^ v {x,y) is its 4D part. Note that the {/i5} and {55} equations, ([7]) and ([8]), 
should be satisfied in the bulk, y ^ , as well as on the brane, y = 0. 

Let us turn to cosmological solutions. The metric is parametrized as follows: 

ds 2 = -P 2 (t, y)dt 2 + Q 2 (t, y)-f ij dx i dx j + £ 2 (t, y)dy 2 . (9) 
The self-accelerated solution reads [2]: 

P{t,y) = l + \y\ JL- , Q{t,y) = a(t) + \y\Va 2 + k, E(t j3 /) = 1. (10) 
Vfl + k 

Here we included a nonzero spatial curvature k. With this Ansatz, the Friedman 
equation on the brane follows from ([3]), and can be expressed in terms of the 4D 
Hubble parameter H = a/a. For a simplest case of k = it reads [2] 

H 2 -m c \H\ = 0. (11) 

This admits a dS solution with H = m c . We emphasize that the minus sign between 
the two terms in ffTTj) is guaranteed by the choice of the positive sign in front of the 
terms in ( jTOl that are proportional to \y\. If we were to choose the latter signs to 
be negative, we would have obtained the Friedman equation H 2 + m c \H\ = 0, which 
does not admit the dS solution. 

On the other hand, as was discussed in section 1, it is exactly the positive sign in 
front of the \y\ terms in (|T0|) that gives rise to the negative mass difficulties. Hence, 
the question: can we avoid these problems while still retaining good properties of 
the self-accelerated solution? 

The proposal is to flip the sign in front of the second term on the l.h.s. of 0- 
In other words, we introduce a new set of equations in which ([3]) is replaced by: 

G^u + m c (K^ - g^K) = MG N T^(x) , (12) 

while all the other equations ©,([7]) and ([8]) remain intact. The action functional 
that gives rise to this new set of equations will be discussed in the next section. 
Here, for the purposes of finding a classical solution and perturbations about it, it 
is enough to focus on the equations of motion. 

The metric for the self-accelerated solution of the new system of equations (1T21 
EHED reads: 

P{t,y) = l-\y\ J— , Q(t,y) = a(t)-\y\Va 2 + k, E{t,y) = l, (13) 
v a ~\~ k 

where we have chosen a negative sign in front of the terms proportional to \y\. Let 
us now see how the change of the positive signs in the metric ffTUl) into the negative 
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signs in (|T3j) changes the value of the extrinsic curvature evaluated at y = + . On 
the solution (IT5|) . = 0, N = 1, and = d y g^ u /2. Hence, at y = + the 
components of the extrinsic curvature tensor evaluated on the solution ([TBI) equal 
to minus the corresponding components evaluated on fflOl) . Therefore, substitution 
of ffTUj) into (jHJ would give the same equation as the substitution of (fT3|) into (TT2T) . 
The corresponding Friedmann equation on the empty brane, which now follows from 
( fT2"j) instead of Q, reads: 



k /„„ fc 



H 2 + ^ = m c ^H 2 + ^. (14) 



For A; = this coincides with ( fill and gives the spatially-flat dS solution with 
if = m c (T5]). For general k the solutions are: 

cfe 2 = (1 - H\y\f {-dt 2 + a 2 {t){d X 2 + S 2 ( X )dn 2 )} + dy 2 , (15) 

where H = m c and k = —1,0,1 corresponds to the open, flat and closed spatial 
slicing of 4D dS space, for which Sk(x) = sinhx, x, sinx, respectiveljU. 

The solution ffl5l) should satisfy all the bulk equations ([6] -[8]), since in the bulk 
it is locally equivalent to Minkowski space. We checked by direct substitution that 
(USD solves Eqs. §\-E> too. 

It is also straightforward to introduce matter/radiation on the brane. Following 
[2] we obtain the Friedmann equation 



which should be amended by the conventional conservation equation for the fluid 
of density p and pressure p: p + 3H(p + p) =0. The latter being a result of the 
matter stress-tensor conservation V^T^ = 0, which can be verified, e.g., by taking 
a covariant derivative of both sides of (|T2|) and using ([7]). 



3 Perturbations 

In this section we study the spectrum of small perturbations of the theory (1121 [6] - 
[S]). Following [9] we consider the metric perturbations of the form: 

ds 2 = (^(x, y) + 5g^(x, y)) dx^dx v + 25g^dx^dy + (1 + Sg 55 (x, y)) dy 2 . (17) 

The background metric will be distinguished by the bar g^ v (x,y) = A 2 (y)^ ^(x) , 
where A = (1 — H\y\), and 7^(x) is a metric tensor for 4D dS space-time. 

3 There are two other solutions to (fl4|) . For k = one finds the H — flat solution. For k = — 1 
one finds the Milne solution a(i) = t. 
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Let us consider small perturbations about an empty (i.e., with = 0) self- 
accelerated background (IT5"|) for k = 0. The following expansion of the metric and 
extrinsic curvature tensors will be used: 

9nv = 9^u + Sg^ , = + 5K pv , (18) 

K^ = ng»v, K = g^K pv = An n=%£, (19) 

where 

K = gTKp, = K + 6K, 5K = g^SK^ - ng^Sg^ . (20) 

It is straightforward to check that the perturbations of the off-diagonal equation (J7J) 
satisfy the following relation 

V (SKfu, - n6 9llu ) = V, (g pa SK p(T - ng p °5g p(7 ) , (21) 

where in this section V denotes a 4D covariant derivative for g^ v . 
Furthermore, small perturbations of the {55} equation (jHJ) yield 

6R = Gnig^SK^ - ng^Sg^) , (22) 

while from the source-free junction condition (fi"2l we obtain 

5R\ 0+ = -3m c (g^5K^ - ng^5g^)\ 0+ . (23) 

Since m c = H, equations fl22|) and ff23l are in contradiction with each other unless 
the r.h.s. of ( |23l) is zero. Requiring also a continuity of 5R we find that the r.h.s. 
of (122p should be zero for arbitrary y: 

g^6K^-ng^6g^ = 0. (24) 

Finally, let us look at small perturbations of the bulk {/xz/} equation. As long as 
T^p = we can use the Gaussian normal coordinates and simultaneously choose the 
following "gauge"0 

= , g^8g^ = , 5g 55 = % s = . (25) 
Under these conditions the bulk {/iz/} equation takes the form: 

d 2 y 5g, u + (D 4 - AH 2 ) 8 9lw = , (26) 

while the junction condition reads as follows: 

- i (D 4 - AH 2 - Hd y ) Sg^\ 0+ = . (27) 

4 In general, this is not an acceptable gauge fixing condition when the brane is held straight 
and the metric is coupled to a non-conformal source. In the latter case one needs to introduce the 
brane bending mode [16] . We'll briefly comment on this below. 
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Having these equations specified, we can calculate the spectrum of Kaluza-Klein 
(KK) modes. It is convenient to use the following decomposition: 




By introducing a new coordinate z and a function u 



in 



dz = 



dy 



f = A 1 ' 2 



(29) 



My) 



u 



we turn fl26|) into a Schrodinger-like equation with a boundary condition set by (T2T1) . 
The former reads 



The spectrum of ( [301 [311) can be computed. The condition for the modes to be 
localized on the brane translates into the requirement, f_°° dzu^z) < oo. 

The above system admits a localized zero-mode solution. For m = 0, we obtain 
u = cexp(— 3H\z\/2) (c being an constant). The zero-mode is normalizable. 

Furthermore, one can check that there are no other normalizable modes satisfying 
both (|30|) and (I3TT) . However, there are an infinite number of continuum states with 
masses m? KK > 9H 2 /4. These are plane- wave normalizable in the z coordinate, and 
their wave-functions are suppressed on the brane because of the 4D EH term [T7| (a 
similar spectrum for a scalar was discussed in [18J). 

Let us compare the above-obtained spectrum to that on the selfaccelerated back- 
ground ([1]), found in |6|. The two spectra are similar to each other for m 2 KK > 9H 2 /4, 
as both of them consist of massive KK modes with suppressed couplings on the 
brane. However, the localized modes differ. On the background ([T]) one finds a 
localized massive mode with ml = 2H 2 , while the localized mode on ([2]) is massless, 
m = 0. In both massive and massless cases there is an enhanced "gauge symmetry" 
for the localized modes: in the massive case the symmetry transformation is [T9] 



(where a is a gauge transformation parameter) while for the massless case this is just 
a conventional 4D reparametrization invariance Sh°^ u (x) = V^( u (x) + V u (^(x). The 
symmetry transformation (I32p is what allows one to gauge away (in the quadratic 
approximation and on an empty brane) the dynamics of the brane bending mode, 
converting the latter into a Lagrange multiplier field However, the symmetry 




(30) 



while the boundary condition takes the form: 




(31) 



6h* (x) = (V M V, + H 2 lt , u )a{x) 



(32) 



6 



(132]) is broken once non-conformal sources are switched on. This leads to: (a) 
propagation of the bending mode with potential ghost-like instability; (b) breakdown 
of perturbation theory making the linearized calculations unreliable [TO]. 

Things are different in the present case, since we have a localized massless mode. 
The gauge symmetry for it, 5h° (x) = + Vj/Cam is preserved when the sources 
are switched on. Then, one should either take into account the brane bending mode, 
or perform calculations with other acceptable gauge conditions, such as e.g., the 5D 
harmonic gauge. The expectation is that this won't introduce an additional on-shell 
degree of freedom. This will be discussed in detail elsewhere. 

4 Modified action 

In this section we will discuss how one should modify the DGP action in order to flip 
the sign in front of the extrinsic curvature terms in the junction condition, i.e., to 
replace ([3]) by (1X21) . Naively, this can be done by flipping the sign of the 5D Einstein- 
Hilbert term in the DGP action. However, this would lead to a bulk graviton with 
a wrong-sign kinetic term. This is not a good course of action. 

Instead, we will introduce an additional 5D Einstein-Hilbert term on the brane, 
with a small negative coefficient. This is acceptable as long as there is a 4D Einstein- 
Hilbert term on the brane with a large positive coefficient. We will show that the 
effect of the new term is twofold: (I) If allows to flip the desired sign in the junction 
condition; (II) It relaxes the constraint on the bulk fundamental scale M*; the latter 
can now take an arbitrary value below Mp\, as long as the coefficient of the brane 
localized 5D Einstein-Hilbert term is tuned appropriately. 

We will discuss these properties first on a simple scalar field example. Then we 
elevate the construction to gravity. 

4.1 Scalar example 

Let us start with the scalar analog of the DGP action [TJ 

/If 2 f M 3 f 

_ _P1 / d^d^ix, 0)) 2 - -f- J d A xdy{d A cj>{x, y)f , (33) 

where the dimensionless scalar field is to mimic a graviton. We impose Z2 sym- 
metry across the y = boundary, and add to the above action the coupling of to a 
source J, also localized on the brane, J d 4 xdy5(y)J(f>. To get the junction condition, 
we integrate the equation of motion obtained from (1331 w.r.t. y, from 0~ to + . The 
resulting equation written at y = + reads: 

- <9 2 0|, =o - m c d y ^\ y=Q+ = J I Ml, . (34) 

Here we introduced the definition m c = 2M^/Mp 1; and for comparison with gravity, 
flipped the overall sign of the equation. 
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Our goal is to modify the action (1331) so that the resulting junction condition 
has an opposite sign in front of the second term on the l.h.s of (1341) . 

For this we introduce an additional term into the action (1331) . This is just a 
5D kinetic term peaked on the brane. To make things tractable, we smear the 
brane, that is, instead of the Dirac function 5(y), we use its regularized version 
5(y) — > 5(y) = 7r _1 e/(y 2 + e 2 ), with e — > 0. The term that we'll be adding to fl33|) 
then reads: 

M 2 f 

— J dxdy5(y)(d A ^x,y)) 2 . (35) 

With this term included the variation of the action 5S = with the appropriate 
boundary conditions gives: 

- (Mp[ - M 2 )5(y)dl<P - Mldl<j> - d y {{Ml - M 2 5(y))d y <p) = J5(y) . (36) 

Next we take the integral of both sides of this equation w.r.t. y from — e to +e, and 
then turn to the limit 

M -> 0, e -»• 0, M 2 /e ~ M 2 5(0) = M 3 > Ml , (37) 

where we keep M fixed, and its value bigger than M*. The resulting equation reads: 



- d 2 J + 2{M M 2 M * ] dy4> = J I Ml, . (3J 



Finally, introducing 



m ^AM^M!l, ( 39) 



where the positive numerical value of m c will be tuned to the Hubble scale today 
m c ~ Hq ~ 10~ 42 GeV, we get the desired junction conditions 

- dl<P + m c d y <P = J /Ml, . (40) 

The second term on the l.h.s. of ( T4T)]) has a sign opposite to the one of the analogous 
term in (|34l) . This accomplishes our goal0. 



5 The Minkowski space propagator following from (j40|) . would have a resonance-like pole on the 
first Riemann sheet, signaling very slow tachyonic instability [20) . The gravitational counterpart 
of this theory, also admits Minkowski space as a solution, in addition to the dS solution we are 
interested in. Perturbations about the former will have a similar pole on the first Riemann sheet, 
showing instability of the Minkowski solution. This instability should be more severe then in the 
scalar case, as the tachyon-like pole here would appear for a tensor state (i.e., the longitudinal 
scalar mode would acquire a ghost-like kinetic term). This is encouraging, since our eventual goal 
is to have only the stable dS background. In the latter case the pole structure of the propagator 
is modified due to the background. The kinetic term of the longitudinal mode should get cured by 
the dS curvature effects. 
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Two comments are in order here. First, the wrong-sign kinetic term (135|) that 
is peaked only on the brane is dominated by the large positive 4D kinetic term in 
( 13"3"j) . proportional to M|j > M 2 . Second, the number of adjusted parameters here 
is the same as in DGP: In ( 133]) one should tune the value of M* such that the ratio 
2M^/Mpj is of order Hq. On the other hand, in the action with the additional term 
([35]) the value of M# < Mpi can be arbitrary, however, for a given value of M* one 
should tune the value of M so that (I3"9"]) is of order H . In the case of gravity, to 
which we'll turn in the next subsection, this has an additional benefit of relaxing 
the constraint on the bulk scale of gravity, M*, which now can be an arbitrary scale 
somewhat lower than Mp\. 



4.2 Gravity 

The DGP gravitational action is [1] 

M 2 f M 3 f 

S = ~t J d " x ^ R ^ j d * xd yV9~sK(g 5 ), (41) 

where g^ u = g^ u (x,0), and g$ refers to the full 5D metric; R and 1Z are the 
four- dimensional and five- dimensional Ricci scalars respectively, and M*, as be- 
fore, stands for the fundamental gravitational scale of the bulk theory The brane 
is located at y = and Z 2 symmetry across the brane is imposed. The boundary 
Gibbons-Hawking term should be taken into account to warrant the correct Einstein 
equations in the bulk. The matter fields, that are also omitted here for simplicity, 
are assumed to be localized on the brane. 

The above action gives rise to the equations of motion ( I3lf6] l8l). As we discussed 
previously, we'd like our equations to be ( fl2l EHS) instead. To reach this goal, we 
proceed as in the scalar case discussed in the previous subsection. 

We smooth out the brane by replacing 5{y) — > S(y), and add the following 
boundary (worldvolume) term to the action: 



m 2 r - m 2 r 

— J d i xdy5{y)^gl ) n{g,) = -— J d 4 x dy 5(y)y/gN (R + K 2 — K 2 U ) ,(42) 



where the r.h.s. of (142H is obtained by using the standard ADM decomposition. The 
total action in the ADM formalism reads: 



S mod = J d*>Xy/giR{gi) + d*x dy ^N(R + K 2 - K\ 

M 2 f 

__y d 4 xdy5(y)^N(R + K 2 -Kl u ) , (43) 

where it is implied that the 4D EH term is also smeared over the same scale as the 
5D term[§. The equations of motion are straightforward to derive from (143j) . The 



3 For the regularization of 4D and 5D EH terms, see, [21] and [22], respectively. 
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{/xz/} equation in the bulk, and {/^5} and {55} equations read as follows: 

G<g = for |y|> c> (44) 
{Ml - M 2 5{y)) (WK,„ - V V K) = , (45) 
{Ml - M 2 5{y)) (R-K 2 + K^K^) = . (46) 

As in the scalar case, we will be looking at this theory in the limit (1371) . The above 
equations reduce to ([ME]) El 

The Israel condition across the brane gets modified because of the new term in 
fj43|) . In the limit fl37j) the junction condition reads: 

2(M 3 - M 3 ) 

GW + M2 ~ 9^K) = 8nG N T^ (x) . (47) 

If M = 0, as in (jHJ), we get back the result © with m c = 2Ml/M^. However, M 
does not have to be zero. For an arbitrary value of M* we tune the value of M so 
that the crossover scale (139]) . which appears in (|47j), is adjusted to the value of the 
present-day Hubble scale m c ~ H ~ 10~ 42 GeV. Hence, ff4T|) recovers (fT2l . 
A few important comments: 

(I) The number of adjusted parameters here is the same as in DGP. In ( |4"TT) 
one should tune the value of M*, so that the ratio 2M 3 /M| 1 is of order H . This 
constrains ~ 100 MeV. On the other hand, in (143]) the value of M* is arbitrary. 
In order to get the right crossover scale for a given value of M* , one should tune the 
value of M so that (139]) is of order H . This has an additional benefit that the bulk 
scale of gravity, M*, can take an arbitrary value M* < Mp\. 

(II) The limit (1371) is a mathematical simplification. In reality the brane will have 
a width e (in simplest cases e ~ M~ x ). If a nonzero brane width effects are kept 
track of, then the junction condition f|4TI) will contain additional terms proportional 
to M 2 . For instance, among other terms, there'll be ones proportional to 

M 2 { 1 \ 

W (G^ - -g, v {K 2 - Kp - 2{K;K av - KK^J . (48) 

Such terms can be neglected as long as M < Mpi. On the other hand, if M ^ Mpi 
the above terms can introduce modifications in the Friedman equation (jllj) . and/or 
need to rescale the overall Newton's constant. 

(III) In the above construction the function F(y) = {Ml — M 2 5{y)) asymptotes 
to Ml away from the brane, \y\ 3> e, while it approaches the negative value, (M 3 — 
M 3 ) < when y — > 0. Being a continuous function, F{y) should pass through zero. 
On the other hand, F{y) is a function that determines the coefficient in front of the 
5D EH term in the action ([4*3]) . This may rise some concerns. Within the brane, i.e., 
for \y\ ^ e the 5D EH term flips the sign, and the concern may be that this could 

7 Eqs. and also flip the overall sign on the brane for our choice of M . This is important 
for understanding the perturbations on the self-accelerated background. 
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lead to some negative-sign kinetic terms. Precisely when the function F(y) flips its 
sign, a large positive-sign 4D EH term (j43p becomes peaked on the brane. The latter 
dominates, as long as Mpj > M 2 . Moreover, the kinetic term of the longitudinal 
modes, which cannot be helped by the 4D EH term, would become normal due to 
the dS background effects, as the calculations of the previous section suggest. 

A related concern may arise in the region where F(y) approaches zero. There, 
the additional helicity states of graviton could become infinitely strongly coupled 
in a naive perturbation theory. However, as in the case of massive gravity, and the 
DGP model, one should expect that a smooth transition to the F{y) — > limit is 
recovered in the full nonlinear theory [23j H] . By now there are a few examples and 
arguments supporting the above expectation, see, e.g., [2U [25J [26] , [TTI [T2] . 

5 Discussions and outlook 

We discussed a model defined by the equations (|T2| MM), or by the action (|43j) . 
and show that it has the self-accelerated solution with decreasing metric in the bulk 
(j2J). The spectrum of small perturbations on the background (j2J) is different from 
the spectrum on ([T]). The lowest mode is a localized massless spin-2 state. The 
low energy theory of this mode exhibits 4D general covariance. Most importantly, 
however, there is no massive localized spin-2 state, the one that causes difficulties 
for the background (CQ) when non-conformal sources are switched on it [9j. All the 
massive states in the present case appear above the threshold, m 2 > 9if 2 /4. 

Nonlinear dynamics of this model is expected to have features similar to DGP, as 
well as different ones. The differences may arise because of the zero mode. On the 
other hand, the KK gravitons form a massive (though light m ~ H) resonance- like 
state which should participate, along with the zero mode, in a one-graviton exchange 
amplitude. This would lead to the vDVZ discontinuity [27] in the linearized theory. 
However, the continuity should be expected to be recovered in the full nonlinear 
theory [4], with small, but potentially measurable deviations from 4D GR [28, 29J 
(see also, [IS EE]). 

Is the obtained background ([2]) acceptable from the point of view of observational 
cosmology? The time evolution in ([2]) is equivalent to that of ([T]). The latter has 
been studied in detail [3], [21]- [10] • This evolution is rather restrictive. For instance, 
the equation of state parameter today, wq = w(z = 0), gets related to the matter 
density Qm, as wo = —1/(1 + £Im) (this is for the k = case [3j). Fitting the 
supernova data with such an expansion favors a low value for Qm, which is in 
tension with the determination of the matter density from CMB observables [40J. 
It has been suggested in [10] to use the open (k = —1) universe to obtain a better 
fit. Forthcoming measurements of the low z supernovae (which are important as 

8 One could consider adding to the DGP action (|4Tj) only the M 2 F(y)(K 2 — K 2 „) term. However, 
this would also change the {55} equation in addition to modifying the junction condition. Then, 
one would need to study the new background solutions. 
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the calibration points) may affect the above considerations. Moreover, the issues of 
determination of the matter content of the Universe may be influenced by the non- 
perturbative phenomena leading to the mass screening [HJ [T2] . The questions of 
how much screening (or anti-screening ) is there, will depend on a particular model 
at hand. These issues have not been understood yet in the present context. 

Finally, an additional benefit of the new term in (H3j) is that it allows to relax the 
constraint on the bulk gravity scale. The latter can take an arbitrary value below 
Mpi. This opens a window for a possible string theory realization of this model, or 
its D > 5 counterparts [UJI12] (f° r earlier proposals see [131 144"] ) . 

With the new term included, the analog of the conventional (Minkowski) solution 
of DGP could also be studied. In the latter case, the function F(y) should not flip 
its sign as it interpolates from its bulk value, M%, down to lower positive value on the 
brane Mf — M 3 (F(y) could emerge as a profile of some 5D scalar field which couples 
to 5D EH term). It would be interesting to study the nonperturbative solutions, 
similar to those of [11], |12j . in this context. 
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